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SUMMARY 


A study  of  the  motion  of  a personal-owner  or  liaison  airplane 
through  the  stall  and  into  the  incipient  spin  has  been  made  by  cata- 
pulting a dynamic  model  into  still  air.  The  results  showed  that  the 
model  became  unstalled  in  inverted  flight  after  its  initial  rotation  of 
substantially  l80°  in  roll.  Since  the  rates  of  rotation  were  low,  it 
may  be  inferred  that  termination  of  the  incipient  spin  would  be  most 
readily  obtained  by  proper  control  movements  at  this  point. 


INTRODUCTION 


Statistics  indicate  that  a large  percentage  of  all  fatal  private 
flying  accidents  in  the  past  have  occurred  because  of  pilot  error;  in 
more  than  half  of  these  accidents,  the  airplane  stall  has  been  involved 
(ref.  l) . Oftentimes,  accidents  of  this  nature  are  referred  to  as 
"stall-spin"  accidents  and  they  occur  primarily  in  the  incipient  phase 
of  a spin,  that  is,  in  that  portion  of  the  motion  immediately  following 
the  stall  and  before  a spin  can  fully  develop.  In  addition,  spin  demon- 
strations of  some  airplanes  for  the  armed  services  require  recoveries 
from  only  two-turn  spins  which  also  may  not  be  fully  developed  motions. 
Data  on  the  incipient  spin  are  not  available  although  much  effort  has 
been  expended  in  the  Langley  20-foot  free-spinning  tunnel  toward  an 
understanding  of  the  developed  spin.  The  problems  associated  with  the 
incipient  phase  of  the  spin  may  be  quite  different  from  those  associated 
with  the  fully  developed  spin  and  recovery  therefrom. 

A technique  is  being  developed  with  free  dynamic  models,  similar 
to  those  used  in  spin  testing,  to  study  the  incipient  spin.  The  model 
is  catapulted  into  still  air  in  such  a manner  that  stalled  flight  occurs. 
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and  the  ensuing  motion  at  and  "beyond  the  stall  is  studied  by  use  of 
motion-picture  records  from  which  the  angles  of  attack  and  sideslip  and 
the  linear  and  angular  velocities  are  determined.  Scale  effect  could 
greatly  influence  the  motion  at  and  Just  beyond  the  stall  although  it 
has  been  shown  to  be  of  only  secondary  importance  in  developed  spins. 

The  techniq^ue  described,  therefore,  may  be  used  primarily  to  isolate  the 
effects  of  various  parameters  rather  than  to  Indicate  specific  incipient- 
spin  characteristics.  The  initial  investigation  was  performed  to  deter- 
mine the  motion  in  the  incipient  spin  of  a model  of  a low-wing,  single- 
vertical-tail  airplane  typical  of  a present-day  four-place  personal-owner 
or  liaison  airplane  and  to  determine  the  extent  to  which  the  motion  may 
repeat  Itself.  The  results  of  this  investigation  are  reported  herein. 


SYMBOLS 


The  analysis  of  the  motion  was  made  with  respect  to  the  body  system 
of  axes . A diagram  of  these  axes  showing  various  angles  pertinent  to 
the  analysis  is  given  in  figure  1. 


x,y,z 


eg 


earth  axes,  a system  of  mutually  perpendicular  axes 
through  airplane  center  of  gravity  always  per- 
pendicular and  parallel  to  earth's  surface  and  of 
constant  direction 

body  axes 

distance  along  X earth  axis  from  reference  point  to 

right  wing  tip,  tail,  and  center  of  gravity,  respectively 

distance  along  Y earth  axis  from  reference  point  to 

right  wing  tip,  tail,  and  center  of  gravity,  respectively 

distance  along  Z earth  axis  from  reference  point  to 

right  wing  tip,  tail,  and  center  of  gravity,  respectively 

dists^nce  of  right  wing  tip  from  center  of  gravity  along 
X,  Y,  and  Z earth  axes,  respectively 

distance  of  tail  from  center  of  gravity  along  X,  Y, 
and  Z earth  axes,  respectively 

span  of  model,  ft 

tail  length  (distance  between  center  of  gravity  and  inter- 
section of  Xb  with  trailing  edge  of  rudder) , ft 
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c 


mean  aerodynamic  chord,  ft 


x/c  ratio  of  distance  of  center  of  gravity  rearward  of  leading 

edge  of  mean  aerodynamic  chord  to  mean  aerodynamic  chord 

z/c  ratio  of  distance  between  center  of  gravity  and  fuselage 

reference  line  to  mean  aerodynamic  chord  (positive  when 
center  of  gravity  is  below  fuselage  reference  line) 


m 


mass  of  airplane,  slugs 


moments  of  inertia  about  Xb>  Yb,  and  Ztq  axes , 

p 

respectively,  slug-feet'" 


mb  2 


inertia  yawing-moment  parameter 


inertia  rolling-moment  parameter 


\ ~ % 
mb^ 

P 

V 

u,v,w 

a 

P 

h 

P 


inertia  pitching-moment  parameter 
air  density,  slugs/cu  ft 

resultant  velocity  of  center  of  gravity,  ft/sec 

components  of  velocity  V along  X,  Y,  and  Z earth 
axes,  respectively,  ft/sec 

components  of  velocity  V along  Xb,  Yb>  and  Zg  axes, 
respectively,  ft/ sec 

Wb 

angle  of  attack,  tan”-*-  — , deg 

^B 

-1 

angle  of  sideslip,  sin”-*-  deg 
altitude  loss,  ft 

rolling  angular  velocity  about  Xb,  radlans/sec 


k 
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q 

pitching  angular  velocity  about  Yg,  radians/sec 

r 

yawing  angular  velocity  about  Zg,  radians/ sec 

« 

OA  and  OC 

projections  of  Xg  and  Yg,  respectively,  on 
earth  X,Y  plane  (fig.  l) 

OB 

position  Xg  would  take  if  airplane  were  rotated  about 
Yg  until  Xg  lay  in  earth  X,Y  plane  (fig.  l) 

OD 

position  Yg  would  take  if  airplane  were  rotated  about 
Xg  until  Yg  lay  in  earth  X,Y  plane  (fig.  l) 

OE 

position  Zg  wo\ild  take  if  airplane  were  rotated  about 
Yg  until  Xb  lay  in  earth  X,Y  plane  (fig.  l) 

OF 

position  Zg  would  take  if  airplane  were  rotated  about 
Xb  until  Yb  lay  in  earth  X,Y  plane  (fig.  l) 

9 

angle  that  Xb  axis  makes  with  earth  X,Y  plane, 
positive  when  nose  is  above  plane  (Euler's  angle  0) 

0 

angle  between  Yb  and  line  OD  (fig.  l)  (Euler's 
angle  0) 

♦ 

angle  between  X earth  axis  and  QA  measured  as  shown 
in  figure  1 from  0°  to  360°  (E-uler's  angle  y) 

t 

rate  of  change  of  with  time,  radians/sec 

rate  of  change  of  0 with  time,  radians/ sec 

0 

rate  of  change  of  0 with  time,  radians/sec 

V 

angle  between  Y earth  axis  and  OC  measured  as  shown 
in  figure  1 from  0°  to  3^0° 

0’ 

angle  that  Yb  axis  makes  with  earth  X,Y  plane  posi- 
tive when  right  wing  is  below  plane 

0' 

angle  between  Xg  and  line  OB  (fig.  l) 

f 

angle  in  X,Y  plane  between  QA  and  OB  or  between  OC 
and  OD  (fig.  l) 

* 
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APPARATUS  AM)  METHODS 
Model 


The  model  used  in  this  investigation  was  a configuration  similar 
with  respect  to  mass  and  dimensional  characteristics  to  one  investigated 
in  reference  2 for  fully  developed  spin  and  recovery  characteristics. 

Based  on  the  model  size  and  the  average  dimensions  obtained  for  a large 
number  of  personal-owner  airplane  designs,  the  model  may  be  considered 
to  be  a 1/12.4  scale  of  a personal-owner  or  liaison  airplane.  The  weight 
and  moments  of  inertia  of  the  model  represent  those  of  an  airplane  typical 
of  this  type.  A three-view  drawing  of  the  model  is  shown  in  figure  2 and 
the  mass  and  dimensional  characteristics  of  a corresponding  full-scale 
airplane  are  tabulated  in  table  I. 

The  model  was  constructed  principally  of  balsa  and  reinforced  with 
spruce  and  cedar.  A photograph  of  the  model  is  shown  in  figure  3*  In 
accordance  with  spin- tunnel  practice,  the  propeller  and  landing  gear  were 
not  simulated  on  the  model.  The  model  wing  had  an  NACA  U3OI2  airfoil 
section  which  at  low  Reynolds  number  has  a shai^p  break  in  the  lift  curve 
at  the  stall.  (See  ref.  3.) 


Apparatus 


The  tests  were  performed  inside  a building  approximately  JO  feet 
square  and  60  feet  high.  The  launching  apparatus,  shown  in  figure  4, 
was  located  55  feet  above  the  floor  near  one  wall  of  the  building  and 
consisted  primarily  of  an  elastic  chord  which  catapulted  the  model  from 
its  launching  platform  along  a short  track . The  launching  platform 
could  be  adjusted  to  set  the  model  attitude  and  thus  the  launching  angle 
of  attack.  An  electronic  timer  was  used  to  measure  the  model  velocity 
Just  after  it  was  la\uiched.  Fixed  motion-picture  cameras  were  located 
along  the  walls  parallel  and  perpendicvilar  to  the  initial  launching 
direction  to  record  the  model  motion.  The  cameras  had  overlapping  fields 
so  that  the  model  was  in  view  of  at  least  two  of  the  cameras,  located  at 
right  angles  to  one  another,  throughout  its  flight.  The  motion-plctxire 
records  were  synchronized  by  use  of  a common  timing  light  reflected  into 
the  camera  fields.  A large  net  for  model  retrieving  purposes  was  hving 
from  the  wall  opposite  the  launching  appeuratus. 


Testing  Technique 

The  model  was  catapulted  from  the  launching  apparatus  at  a speed 
slightly  in  excess  of  the  stalling  speed  and  at  em  angle  of  attack  Just 
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"below  the  stall  angle  of  attack.  The  elevator  control  of  the  model  was 
set  so  as  to  pitch  the  model  from  the  launching  angle  of  attack  up  through 
the  stall  and  the  rudder  was  set  to  initiate  a yawing  motion  and  thus  to 
precipitate  a roll-off  after  stalling. 

As  previously  mentioned,  photographic  time  histories  of  the  model 
motion  were  made  with  a system  of  motion-picture  cameras  arranged  so  as 
to  have  the  model  in  the  field  of  two  cameras  with  perpendicular  focal 
axes  at  all  times . When  at  an  instant  of  time  the  model  is  in  the  field 
of  two  cameras  located  at  right  angles  to  one  another,  the  space  coordi- 
nates of  the  model  can  he  determined.  By  this  means,  the  model  could  be 
located  at  any  instant  from  the  time  of  catapult  until  the  time  the  model 
struck  the  retrieving  net.  From  a knowledge  of  the  space  coordinates 
and  their  instantaneous  time  values,  the  velocities  of  the  model  were 
obtained  along  the  earth  axes . Then,  the  attitude  of  the  model  at  any 
instant  being  known,  angle  of  attack,  angle  of  sideslip,  and  linear  and 
angular  velocities  were  calculated.  The  derivation  of  equations  for  these 
calculations  is  presented  in  the  appendix. 


Test  Conditions  and  Accuracy 

For  the  tests,  the  model  was  ballasted  with  lead  weights  to  repre- 
sent a corresponding,  airplane  at  an  altitude  of  5,000  feet  (p  = 0.0020i|-9 
slug/cu  ft) . 

For  this  investigation,  the  controls  of  the  model  were  set  so  as  to 
cause  the  model  to  stall  and  enter  a spin  soon  after  launching.  The 
settings  of  the  controls  were  as  follows: 


Elevator  full  up,  deg 30 

Ailerons,  deg 0 


Rudder  fully  deflected  (in  direction  desired  for  roll-off),  deg  . . . 30 

The  stalling  angle  and  stalling  speed  of  the  model  were  determined 
from  force  tests  of  the  model  conducted  on  a strain-gage  balance  in  the 
Langley  20-foot  free-spinning  tunnel.  The  Reynolds  number  of  the  present 
investigation  was  approximately  1 x 10^  based  on  the  mean  aerodynamic 
chord  of  the  model. 

Errors  in  the  measured  motion  were  magnified  because  graphical  dif- 
ferentiation was  required  in  reducing  the  data.  For  the  c\irrent  tests, 
the  acciiracy  was  estimated  to  be  within  the  following  limits: 
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a,  deg ±5 

P,  deg ±4 

u,  V,  w (full  scale),  ft/sec -1-3 

uj,  Vg,  wgCfiill  scale),  ft/sec ±4.9 

p,  q,  r (full  scale),  radians/sec ±0.15 

Altitude  loss,  (full  scale),  ft_ ±1 

Weight,  percent  ±1 

Center-of -gravity  position,  percent  c ±1 

Moments  of  inertia,  percent  ±5 

Control  settings,  deg  ±1 


RESULTS  AM)  DISCUSSION 


The  results  of  the  tests  are  presented  In  figures  5 and  6 . The 
model  test  data  have  "been  converted  to  full-scale  values  for  a corre- 
sponding airplane  at  a test  altitude  of  5,000  feet  "based  on  the  premise 
that  the  model  scale  was  l/l2.4.  Film  strips  showing  the  motion  of  the 
model  as  it  stalls  and  enters  a spin  are  presented  in  figure  J. 

The  results  as  the  model  rolled  to  the  right  and  entered  an  incipient 
spin  are  presented  in  figure  5 for  an  initial  and  a repeat  test.  Varia- 
tions with  time  of  the  angle  of  attack,  angle  of  sideslip,  altitude  loss, 
and  the  rates  of  roll,  pitch,  and  yaw  are  given  in  figure  5.  Sketches  of 
the  approximate  attitude  of  the  model  at  various  phases  of  the  motion  are 
also  shown  in  figure  5.  As  was  previously  noted,  the  elevator  was  full 
up  (30°)  and  the  rudder  was  full  right  (30°)  for  these  tests.  Soon  after 
the  model  stalled,  it  again  "became  unstalled  (time,  approx.  1.75  sec), 
went  inverted  for  a part  of  its  motion,  and  assumed  a very  low  angle  of 
attack  momentarily  before  again  stalling  and  continuing  toward  a developed 
spin  motion,  such  as  reported  in  reference  2.  After  the  initial  stall  and 
immediately  after  the  model  "became  unstalled,  the  rates  of  yaw  and  pitch 
were  relatively  small  and  the  rate  of  roll  "began  to  decrease.  There  also 
was  little  loss  of  altitude  to  this  time  and  the  results  indicated  that, 
even  thovigh  the  airplane  may  he  inverted,  a time  soon  after  the  roll-off 
has  started  appears  to  be  a desirable  time  to  attenpt  to  terminate  the 
motion.  This  conclusion  is  in  agreement  with  pilots ' experience  that 
recoveries  attempted  soon  after  an  airplane  stalls  and  rolls  off  are 
more  readily  obtainable  than  those  attempted  from  developed  spins. 

The  motion  obtained  may  be  qxialitatively  analyzed  as  follows:  After 

being  launched,  the  model  immediately  pitched  up  euid  yawed  right  because 
of  the  initial  settings  of  the  elevator  and  rudder,  respectively.  As  the 
motion  progressed,  the  sideslip  angle  p became  negative  becavise  of  the 
yawing  motion;  this  sideslip  euagle  tended  to  cause  the  model  to  roll  off 
as  a resvilt  of  the  dihedral  effect.  As  was  previovtsly  noted,  the  wing 
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had  an  NACA  45012  airfoil  section  for  which  the  lift  curve  has  a sharp 
break  at  the  stall  (ref.  5)»  this  effect  undoubtedly  augmented  the  model 
roll-off  because  of  a loss  in  damping  in  roll  beyond  the  stall.  When 
the  model  again  became  unstalled,  the  damping  in  roll  led  to  decreased 
rate  of  roll  (p)  until  the  stall  reoccurred  with  the  resulting  increase 
in  rate  of  roll.  Changes  in  sideslip  euigle  occurred  as  the  model  rolled 
and  Increments  of  positive  and  negative  velocity  components  along  the 
Yb  axis  were  obtained  as  the  right  or  left  wing,  respectively,  was 
Inclined  downward. 

One  of  the  main  objectives  of  the  investigation  was  to  ascertain 
whether  motions  through  the  stall  and  into  the  incipient  spin  could  be 
reliably  indicated  by  means  of  free  dynamic  models.  The  repeat  test, 
the  results  of  which  are  shown  in  figure  5,  was  made  with  the  test  con- 
ditions nearly  the  same,  insofar  as  the  available  equipment  allowed, 
as  those  of  the  original  or  initial  test  in  order  to  determine  whether 
the  motion  would  repeat  itself.  Although  there  were  some  differences 
for  the  repeat  test  in  that  the  model  had  a more  positive  pitching 
velocity  and  a somewhat  different  sideslip  angle  than  for  the  initial 
test,  the  results  in  general  are  in  good  agreement.  Quantitatively, 
the  results  are  generally  close  to  or  within  the  estimated  limits  of 
accuracy  and  indicate  that  the  model  probably  repeats  its  motions  for 
similar  launching  conditions.  Because  of  the  unsteady  nature  of  the 
air  flow  which  exists  at  and  beyond  the  stall,  it  has  been  tho\ight  that 
airplanes  may  not  repeat  motions  in  the  stalled- flight  regime.  The 
results  from  the  repeat  test  and  experience  in  spin  and  spin-recovery 
work,  however,  indicate  that  motions  at  and  beyond  the  stall  may  gen- 
erally repeat  themselves  except  for  critical  cases  such  as  have  been 
encountered  in  studies  of  developed  spins  when  the  spin  and  recovery 
change  vinaccountably . 

The  motions  discussed  are  for  power-off  stalls  emd  the  use  of  power 
may  appreciably  alter  the  motion  or  the  symmetry  from  that  presented 
herein.  The  resvilts  should  not  be  Interpreted  as  representing  the 
motions  of  an  airplane  corresponding  precisely  to  the  model,  because  of 
scale  effects,  but  can  be  considered  indicative  of  the  trends  in  the 
motion  of  a typical  personal-owner  or  liaison  airplane.  Changes  in 
configuration,  wing  loading,  and  mass  distribution  also  may  alter  the 
results  appreciably. 

In  order  to  evaluate  the  symmetry  of  the  model  and  to  ascertain 
that  the  motion  to  the  right  previously  discussed  was  not  Influenced 
primarily  by  model  asymmetry,  tests  were  conducted  to  determine  the 
motion  with  the  rudder  set  full  left  (30°) . The  results  are  presented 
in  figure  6.  As  before,  the  variation  of  angle  of  attack,  angle  of 
sideslip,  and  angular  velocities  with  time  are  presented  as  the  model 
stalled  and  rolled  to  the  left  prior  to  entering  a spin  to  the  left. 
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A con^jarison  of  figure  6 with  figure  5 indicates  that  the  model,  which 
appeared  to  be  dimensionally  symmetrical,  showed  generally  similar  motions 
when  the  rudder  was  set  either  to  the  right  or  to  the  left. 


CONCLUDING  REMARKS 


An  investigation  has  been  conducted  to  study  the  motion  of  a free 
dynamic  model  of  a personal-owner  or  liaison  airplane  through  the  stall 
and  into  the  incipient  spin.  This  dynamic-model  technique  does  not  lend 
itself  to  specific  or  development  studies  of  a given  design  because  of 
probable  scale  effects  existing  at  the  stall;  however,  the  method  is 
applicable  for  Indicating  the  trends  of  the  motion  through  the  stall  and 
into  the  incipient  spin  and  possibly  the  effect  of  certain  parameters  on 
this  motion.  The  results  showed  that  the  model  became  unstalled  in 
inverted  flight  after  its  initial  rotation  of  substantially  l80°  in  roll. 
Since  the  rates  of  rotation  were  low,  it  may  be  inferred  that  termination 
of  the  incipient  spin  would  be  most  readily  obtained,  by  proper  control 
movements  at  this  point.  The  results  of  the  model  flights  show  that 
the  motion  is  generally  repeatable  and  that,  power  effects  being  excluded, 
a geometrically  symmetrical  model  should  have  fairly  symmetrical  motions 
to  the  right  and  to  the  left. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  January  28,  1953* 
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APPENDIX 

DERIVATION  OF  EQUATIONS  FOR  EVALUATION  OF 
PARAMETERS  IN  INCIPIENT-SPIN  MOTION 


The  determination  of  angle  of  attack,  angle  of  sideslip,  and  the 
angular  velocities  of  the  model  in  flight  are  based  upon  calculating 
its  component  velocities  along  the  body  axes.  The  determination  of 
these  velocities  is  accon^lished  by  differentiating  the  respective 
earth-axis  space  coordinates  of  the  model  against  time  to  obtain  the 
earth-axes  velocities  of  the  model  and  then  resolving  by  trigonometric 
means  these  velocities  to  a body-axes  system.  The  following  process 
was  used. 

If  X,  Y,  and  Z are  the  eeurth  axes,  and  Xp.  Yg,  and  Z,g  are 
the  body  axes,  both  systems  having  their  origin  at  the  center  of  gravity 
of  the  model,  the  relationships  between  earth-axes  velocities  and  body- 
axes  velocities  are  as  follows: 


Ug  = u cos  l_  XXg  + V cos  l_  YXg  + w cos  [_  ZXg  (A1) 


Vg  = u cos  !_  XYg  + V cos  !_  YYg  + w cos  f_  ZYg  (A2) 


Wg  = u cos  l_  XZg  + V cos  j_  YZg  + w cos  j_  ZZg  (A3) 


where  / XXg,  / YXg,  and  so  forth,  are  angles  between  the  earth  and 

body  axes  designated.  These  nine  angles,  which  relate  the  velocities 
along  the  earth  and  body  axes  throiogh  direction  cosines,  therefore  must  be 
found.  The  coordinates  of  the  model  right  wing  and  tail  reference  points 
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(wing  tip  In  Xg,Yg  plane  and  tail  cone  in  plane,  respectively) 

from  the  center  of  gravity  are  determined; 


^ ^cg 


= ^-R  - ^-cg 


K Tp  r^g 


^ ’ ^cg 


yT  = yrj  - yr^g 

2T  = 2r,j,  - Zr^jg 


From  the  foregoing  coordinates,  the  angles  0,  0',  i|r',  0', 

and  0,  as  well  as  ^r”,  can  he  obtained  (see  fig.  l) ; 


0 = sin  — 

Zip 

1 yT 

* = tan"^  — 

Xip 


0-  =sin-^!^ 

b/2 

i|r'  = tan"^ 

yR 

♦”  = t'  - t 


0*  = C08“1(C0S  0 cos  if") 


(Ak) 


0 * COS"l(cOS  0'  cos  f") 
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From  figures  1 and  8 and  from  the  determination  of  the  angles  cal- 
culated in  equations  (a4) . the  trigonometric  functions  necessary  for  the 
solution  of  equations  (A1)  and  (A2)  can  he  obtained  as  follows: 


cos  /XXb  = cos  \|f  cos  e (A5) 

cos  / YXb  = cos(90  - y)cos  9 = sin  \|f  cos  6 (a6) 

cos  / ZXb  = cos (90  + e)  = -sin  0 (AT) 

cos  J_  YYb  = cos  V cos  0'  (a8) 

cos  Zxyb  = cos(90  + \|f')cos  0'  = -sin  \|r ' cos  0'  (A9) 

cos  L zyb  = cos(90  - 0')  = sin  0'  (AlO) 


The  determination  of  the  trigonometric  functions  needed  for  the 
solution  of  eqviation  (A3)  is  somewhat  more  difficult  than  those  for 
equations  (A1)  and  (A2) . These  functions  may  he  found  hy  consideration 
of  the  spherical  triangles  depicted  in  figure  8. 

cos  ZZb  = cos  9 cos  0 = cos  0'  cos  0'  (All) 

Also,  from  figure  8 

cos  Z XZg  = sin  !_  ZZg  cos  j_  XZZb 
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where 


and 


sin  j_  ZZg  = ITT 


cos^e 


cos 


20 


cos  !_  XZ2/Q  = cos  cos  /_  2,gZF  + sin  \lr  sin  / 


or 


cos  j_  XZZg  = 


cos  \|r  sin  Q cos  0 + sin  sin  0 
'1  - cos^e  cos20 


Therefore, 


cos  !_  X2ig  = cos  \lr  sin  9 cos  0 + sin  \|f  sin  0 


(A12) 


Now 


cos  / YZg  = sin  / ZZg  cos  £ YZZb 


where 


or 


cos  /_  YZZg  = sin  cos  £ 2igZF  - cos  if  sin  £ ZgZF 


cos  £ YZZg  = 


sin  y Sin  9 cos  p 


\ 1 - COS^e  cos^0 


Ik 
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Therefore, 


cos  / YZg  = sin  t sin  0 cos  0 - cos  if  sin  0 (AI3) 

(Note  that  angles  XZZg,  YZZg,  and  ZgZF  are  angles  between  the  arcs 

of  the  spherical  triangles  shown  in  figure  8.)  Thus,  with  the  equa- 
tions (as)  to  (AI3),  the  nine  necessary  direction  cosines  are  determined 
and  eqmtions  (Al)  to  (A3)  may  be  solved  for  the  velocities  along  the 
body  axes.  The  velocities  along  the  earth  axes  necessary  for  the  solu- 
tion of  these  equations  are  obtained,  as  was  indicated  previovisly,  by 
taking  the  first  derivative  of  the  space  coordinates  with  respect  to 
time. 


From  a solution  of  equations  (Al)  and  (A3)  wherein  a determination 
of  the  velocities  of  the  center  of  gravity  along  the  X and  Z body 
axes  has  been  made,  the  angle  of  attack,  measured  at  the  center  of  gravity 
and  defined  as  the  angle  between  the  X body  axis  and  the  projection  of 
the  relative  wind  on  the  Xb,Zb  plane,  can  be  calculated  as  follows: 


a 


tan”l 


(AlU) 


If  the  resultant  velocity  of  the  model  is  calculated 


V 


(A15) 


the  of  sideslip  may  be  calcvilated  as  follows: 


P 


(AI6) 


For  a determination  of  the  euigular  velocities  p,  q,  and  r,  refer- 
ence 4 was  used.  The  angular  velocities  are  determined  on  the  basis  of 
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Euler's  angles,  which  are  herein  defined  as  if,  9,  and  0.  The 
tlons  for  the  angular  velocities  therefore  are  as  follows; 


P 


= 0 - ^ sin  9 


q = 9 cos  0 + ijr  sin  0 cos  9 


eq.ua- 


(A17) 


(Al8) 


r = cos  0 cos  9-9  sin  0 


(A19) 
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TABLE  I.-  DltOSIOHAL  ADD  MASS  CHARACTERISTICS  OF  THE 
CORRESPOMDIRa  FULL-SCALE  AIRPLANE 


Airfoil  section  NACA  43012 

Incidence,  deg 3 

Dihedral,  deg 6 

Twist,  deg 0 

Span,  ft 33.63 

Mean  aerodynamic  chord,  c,  ft 4.89 

Leading  edge  of  mean  aerodynamic  chord  rearward  of  leading  edge  of  wing,  ft 0.05 

Taper  ratio  . . . 1.00 

Area,  sq  ft 163.28 

Aspect  ratio  6.93 

Ailerons ; 

Span,  ft 7,15 

Area,  rearward  of  hinge  line,  sq  ft 15,70 

Horizontal-tall  surface: 

Span,  ft 10.25 

Total  area,  sq  ft 26.39 

Elevator  area  rearward  of  hinge  line,  sq  ft 11.02 

Aspect  ratio  3.98 

Incidence,  deg 0 

Dihedral,  deg _ 0 

Distance  from  quarter  chord  of  c to  elevator  hinge  line,  ft  13.73 

Section Modified  NACA  0009 

Vertical-tail  surface: 

Span,  ft 5.32 

Total  area,  sq  ft 12.^ 

Rudder  area  rearward  of  hinge  line,  sq  ft 6.48 

Aspect  ratio  1,26 

Offset,  deg 0 

Dlstsmce  from  quarter  chord  c to  rudder  hinge  line,  ft  l4.l8 

Section Modified  NACA  0009 

lass: 

Weight,  ih 2449 

Wing  loadii>g,  Ih/sq  ft 15.0 

Relative  density: 

Sea  level 5.82 

5,000  ft 6.76 

Center  of  gravity: 

x/c 0.26 

z/c 0,12 


Moments  of  Inertia  (about  center  of  gravity),  slug-ft^ 

^Xb 



^Zb 


Inertia  parameters: 

% ‘ % 

It-  - Izb 


ii».S 


I 
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Figure  1.-  Diagram  of  earth  and  body  axes  through  the  center  of  gravity 
of  the  model  showing  angular  relations.  Positive  directions  of  axes 
and  positive  values  of  angles  are  shown  except  for  9 and  9 ' which 
are  shown  negative. 
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Figure  2,-  Three-view  drawln 


Incidence 


nuwvT 

Nnge  line 


of  model  tested 
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Figure  3.-  Photograph  of  model  tested. 
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(a)  Angle  of  attack,  angle  of  sideslip,  and  edtltvide  loss. 

Figure  5.-  Tine  histories  of  sone  coiqponents  of  model  notion  for  the 
original  and  repeat  right  roll-offs. 
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(b)  Bolling,  pitching,  and  yawing  velocities. 
Figure  5.-  Concluded. 


2 3 4 5 6 

Time,  «ec 


(a)  Angle  of  attack,  angle  of  sideslip,  euid  altltiide  loss. 

Figure  6.-  Time  histories  of  some  conQ>onents  of  model  motion  for  the 

symmetry-check  left  roll-off. 


Pitching  velocity , q ,rodlans/sec  Yawing  velocity  , r , rodions/sec 
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(b)  Rolling,  pitching,  and  yawing  velocities. 
Figure  6.-  Concluded. 
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